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Guided  ion  beam  kinetic  energy  thresholds  in  the  ion-  molecule 
reactions  M+  +  H2  MH+  +  H.  where  M+  =  B+.  A1+,  and  Ga"*"  exceed  by  0.4  to  ca.  5 
eV  the  thermodynamic  energy  requirements  or  theoretically  computed  barrier 
heights  of  these  reactions.  In  addition,  the  formation  of  MD***  occurs  at  a 
significantly  lower  threshold  than  MH+  when  M+  reacts  with  HD.  Moreover,  the 
measured  reaction  cross-  sections  for  production  of  MH'*'  product  ions  are  very 
small  (10*^^  to  10‘2®  cm^).  These  facts  suggest  that  a  "dynamical  bottleneck" 
may  be  operative  in  these  reactions.  In  this  work,  the  eigenvalues  of  the  mass- 
weighted  Hessian  matrix,  which  provide  local  normal-  mode  frequencies,  are 
used  to  identify  locations  on  the  ground-state  MH2'**  potential  energy  surfaces 
where  collisional-  to-  internal  energy  transfer  can  readily  take  place.  In 
particular,  the  potential  energies  at  geometries  where  eigenvalues 
corresponding  to  inter-  fragment  and  to  internal  motions  undergo  avoided 
crossings  are  related  to  the  kinetic  energies  of  apparent  reaction  thresholds. 
This  near-  resonance  energy  transfer  model,  applied  to  M'*’  +  HD  reactions, 
displays  the  experimentally  observed  preference  to  form  MD***  at  lower  collision 
energies  than  MH***  as  well  as  the  fact  that  reaction  thresholds  may  greatly 
exceed  thermodynamic  energy  requirements.  This  model  explains  the  small 
reaction  cross-  sections  in  terms  of  high  energy  content  and  subsequent 
dissociation  of  nascent  MH'*’  (or  MD'*’)  ions.  Although  the  mass-  weighted 
Hessian  matrix  is  used  as  a  tool  in  this  analysis,  the  model  put  forth  here  is  not 
equivalent  to  a  reaction-  path  Hamiltonian  dynamics  approach. 


Uakauiad 


93-02119 


93  ? 


c 


a 


I.  Introduction 


Guided  ion  beam  measurements^  of  cross-  sections  for  reactions  of 
closed-  shell  B+.  A1+,  and  Ga+  ions  with  closed-  shell  'Xg'*’  H2.  Dj,  and  HD  have 

shown  features  that  require  interpretation: 

(iK,rhe  apparent  thresholds  (i.e.,  the  collision  energies  where  product  MH'*'  or 
ions  are  first  formed)  exceed  the  minimum  thermodynamic  energy 
requirements  by  significant  amounts  (e.g..  by  up  to  5  eV  for  Ga***). 

(ii)  In  experiments  with  HD,  MD***  formation  displays  a  lower  energy  threshold 
than  MH'*’. 

(iii)  The  cross-  sections  arc  small  (10**^-  10*20  cm^),  and  are  smallest  for  Ga"** 
and  largest  for  B'*’. 

In  the  present  work,  we  report  findings  that  relate  to  these  experiments  and 
which  allow  an  intepretation  of  much  of  the  dau  in  terms  of  features  of  the  M'*' 
H2  potential  energy  surfaces  in  regions  of  strong  mode  coupling. 

In  particular,  a  mass-  weighted  Hessian  analysis  of  the  local  natural 
frequencies  for  intra-  fragment  and  inter-  fragment  motions  in  regions  of 
strong  repulsive  interactions  shows  that  energy  transfer  may  be  the  rate 
limiting  step  in  these  reactions.  The  collision  energies  needed  to  access 
geometries  where  such  dynamical  resonances  occur  are  correlated  with 
observed  reaction  thresholds.  Moreover,  for  M***  -i-  HD  collisions,  energy 
transfer  that  results  in  formation  is  shown  to  occur  at  lower-  energy  than 

that  producing  MH'*’  which  involves  a  higher-  energy  resonance. 

In  all  cases,  the  electronic  energies  in  such  regions  of  strong  coupling 
approach  (for  B*^)  or  even  exceed  (for  Al'*’  and  Ga'*’)  the  dissociation  energy  of 
H2.  As  a  result,  collisions  that  access  such  regions  produce  MH'*’  or  MD***  with  a 
large  amount  of  internal  vibrational  energy.  In  fact,  these  product  ions  are 
likely  to  possess  enough  internal  energy  to  fragment,  thereby  reducing  the 
MH'*’  (MD'*’)  yield  and  the  measured  reaction  cross-  section  (least  so  for  B'*’  and 
most  so  for  Ga'*’). 
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The  present  use  of  eigenmndes  of  the  mass-  weighted  Hessian  matrix 
differs  from  that  embodied  in  the  so-  called  reaction  path  Hamiltonian^ 
approach.  In  our  model,  the  critical  geometries  do  not  lie  on  the  reaction  path 
and  have  energies  far  in  excess  of  such  a  path  or  of  corresponding  first-  order 
saddle  points  (i.e.,  transition  states).  Our  critical  geometries  relate  more  closely 
to  those  that  are  realized  in  the  experiments'  very  non-  equilibrium  high- 
energy  ion-  molecule  bimolecular  collisions  in  which  the  reagents  possess 
little  internal  energy. 

In  Sec.  II.,  we  describe  the  computational  methods  used  to  compute  the 
potential  energy  surfaces,  gradient  vectors,  and  mass-  weighted  Hessian 
matrices  use  in  this  work.  In  Sec.  III.  we  present  and  discuss  our  potential 
energy  surfaces  and  the  reaction  energetics  they  imply,  and  in  Sec.  IV.  we 
introduce  a  dynamical  model  to  simulate  the  early  stages  of  the  M***  +  H2  (D2  or 
HD)  collisions.  Sec.  V.  describes  our  primary  findings  and  their  relation  to  the 
experimental  data,  and  in  Sec.  VI.  we  summarize. 

II.  Computational  Methods 


A.  Basis  Sets 


For  the  B'*'  +  H2  and  Al*^  H2  calculations,  the  H  atom  basis  was  Dunning's 
augmented  correlation  consistent  (cc)  polarized  valence  triple-zeta  (p-VTZ) 
[5s2pldl  3s2pld]  set  of  functions^.  For  the  B'*'  ion,  the  Dunning  [10s5p2dl  4s3p2d] 
augmented  cc  p-VTZ  basis  set^  was  used,  and  a  total  of  S5  contracted  Gaussian- 
type  basis  functions  resulted  for  BH2'^.  For  the  Al*^  ion,  the  McLean-Chandler 
[12s9pl  6s5p]  basis  sef^  augmented  with  one  3d  polarization  function  (exponent 
0.4)  was  used,  and  total  of  S7  contracted  Gaussian-type  basis  functions  resulted 
for  AIH2'*’.  In  the  case  of  Ga'*',  the  so-  called  SBK  psuedopotential^,  which  treats 
Is,  2s,  and  2p  orbitals  implicitly  and  3s,  3p,  3d,  4s  and  4p  orbitals  explicitly,  was 
used  with  a  (8L,  6dl  4L,  3d)  basis.  For  GaH2'*',  a  6-31 IG**  basis^  was  employed  for 
each  H  atom,  thus  giving  a  total  of  46  explicit  atomic  orbitals. 


B.  Electronic  Configurations  and  Wavefunctions 


>VM  I4S 


In  generating  the  potential  energy  surfaces,  optimal  geometries,  and 
local  harmonic  vibrational  frequencies  reported  here,  the  complete  active 
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space  (CAS)  -based  multiconfigurational  self-consistent  field  (MCSCF)  method 
was  used  to  treat  correlations  among  the  valence  electrons  of  the  MHH'*'  system. 
The  six  valence  orbitals  are  all  those  derived  from  the  metal  ns,  np  and  the  two 
H  Is  orbitals.  The  final  electronic  energies  at  critical  (i.e.,  optimal)  geometries 
were  evaluated  at  the  quadratic  configuration  interaction  including  single, 
double,  and  approximate  triple  excitations  (QCISD(T))  level  to  obtain  more 
quantitative  estimates  of  thermodynamic  data.  In  a  few  situations,  convergence 
difficulties  arose  in  implementing  the  QCISD(T)  calculations,  so  we  resorted  to 
fourth  order  Mpller-  Plesset  perturbation  theory  (MP4)  for  computing  our  most 
accurate  energies. 

As  discussed  in  our  earlier  work  on  no  single  electronic 

configuration  can  describe  even  the  ground  state  of  these  systems  throughout 
C2v  •  or  Cs  reaction  paths.  For  this  reason,  multiconfigurational  methods 

were  required.  In  the  MCSCF  calculations,  the  four  valence  electrons  were 
distributed,  in  all  ways  consistent  with  overall  spatial  and  spin  symmetry, 
among  the  6  valence  orbitals.  This  process  generated  41  electronic 
configurations  of  ^Ai  symmetry  in  the  C2v  point  group  and  65  of  ^A'  symmetry 
in  the  Cs  point  group;  it  yielded  41  configurations  in  the  C„v  group. 

The  above  MCSCF  calculations  on  and  AIH2'''  were  employed,  along 

with  our  Utah  MESSKit^^  analytical  energy  derivative  and  potential  energy 
surface  'walking'  algorithms'^  to  find  and  characterize  (via  geometry  and  local 
harmonic  vibrational  frequencies)  the  local  minima,  transition  states,  and 
reaction  paths  discussed  below.  For  GaH2'''.  we  used  the  GAMESS  program  suite^. 
which  uses  finite-  difference  methods  to  compute  the  Hessian  matrix  from 
analytical  energy  gradients.  The  QCISD(T)  and  MP4  energies  were  computed 
using  the  GAUSSIAN  92  program^. 

III.  Reaction  Energetics 

A.  Potential  Energy  Surfaces 

kI«^2v  Surfaces 

In  Figs.  lA-lC  arc  shown  contour  potential  energy  surfaces  for  the  CAS- 
MCSCF  ground  electronic  states  (which  have  singlet  spin  and  totally  symmetric 
spatial  symmetry)  of  the  three  H2  reactions  considered  here  within  C2v 


symmetry.  The  axes  in  the  graphs  are  R.  the  distance  in  A  from  the  M  nucleus 
to  the  midpoint  of  the  H-H  moiety,  and  r.  the  distance  between  the  two  H  nuclei 
(see  Fig.  2).  The  similarities  among  the  three  surfaces  are  striking,  with  the 
primary  differences  being  results  of  (i)  B‘*‘  being  smaller  than  Al'*'  and  Ga'*',  and 

(ii)  the  H-B-H  bonds  being  stronger  than  the  H-Al-H  bonds  which  arc  a  bit 
stronger  than  the  H-Ga-H  bonds. 

In  each  of  these  surfaces,  four  regions  are  noteworthy: 

(i)  the  asymptotic  region  (R  >  3  A  and  r  near  0.7  A)  where  a  narrow  entrance 
channel  governs  the  approach  of  M'*'  to  H2  and  where  the  energy  variation 
along  the  r  coordinate  is  essentially  that  of  an  isolated  H-H  bond  while  that 
along  R  is  rather  weakly  increasing  as  R  decreases, 

(ii)  the  H-M-H'*’  linear-  ion  region  near  R  =  0  pertaining  to  the  locally  stable 

ion  (for  HAl+H  and  HGa+H  this  ion  is  metastable  with  respect  to  H2  +  A1+  or 

H2  +  Ga'*’;  for  HBH'*'.  the  ion  lies  below  H2  +  B*), 

(iii)  the  "barrier”  connecting  the  entrance  channel  and  the  linear-  ion 
minimum  (the  barrier  regions  are  marked  by  X  in  Figs.  1;  we  refer  to  them  as 
barriers  rather  than  transition  states  because,  as  discussed  later,  they  are 
second-  order  saddle  points  on  these  surfaces), 

(iv)  the  region  of  strong  interaction  where  both  R  and  r  are  relatively  small  as 
a  result  of  which  the  couplings  among  the  internal  modes  are  strong  (see  the 
regions  maiiced  by  Y  in  Figs.  1). 

2.  Collinear  Approach  Surface 

In  Fig.  3  is  shown  a  potential  energy  contour  surface  (as  a  function  of 
the  distance  rMH  between  the  metal  and  the  closest  H  atom  and  r)  for  collinear 
approach  of  B'*’  to  H-H;  the  collinear  surfaces  for  Al'*'  and  Ga'*'  display  similar 
features.  We  found  that  as  rMH  decreases  from  its  asymptotic  value,  the  bending 
vibrational  frequency  at  such  collinear  geometries  is  imaginary  and  its 
magnitude  increases  as  imh  decreases.  Of  course,  as  the  angle  between  the  H-H 
axis  and  the  vector  connecting  M  to  the  center  of  the  H-H  moiety  changes  from 


0®  to  90  ®  (i.c.,  from  collinear  to  C2v  geometry),  the  frequency  corresponding  to 
this  motion  becomes  real,  reflecting  the  stability  of  the  C2v  regions  of  the 
surfaces.  The  negative  curvature  along  the  bending  coordinate  is  caused  by  the 
presence  of  low-  lying  2pn  orbitals  on  B'*'  which,  upon  bending  away  from 

collinear  geometry,  mix  with  and  lower  the  energies  of  occupied  valence 
orbitals  thereby  lowering  the  total  energy. 

The  fact  that  the  potential  surface  becomes  more  and  more  unstable  to 
rotating  the  H-H  bond  axis  away  from  the  M'*’  ion  as  rMH  decreases  means  that 
flux  incident  toward  such  collinear  approaches  will  be  moved,  by  forces 
directed  away  from  linear  geometries,  toward  the  "insertive"  C2v  type 
geometries.  For  this  reason,  it  is  unlikely  that  linear  or  near-  linear 
orientations  play  important  roles  in  the  M***  •fH2-»MH'*'  -t-H  reactions  even 
though  there  is  no  barrier  along  such  paths  in  excess  of  the  thermodynamic 
energy  difference  (calculated  here  to  be  61.  91.  and  94  kcal/mol  for  B***,  Al'*'.  and 
Ga*^.  respectively).  It  is  for  this  reason  that  we  focus  the  majority  of  our  study 
and  analysis  on  the  C2v  (and  near)  pathways  although  these  paths  do 
experience  barriers  in  excess  of  thermodynamic  requirements. 

B.  Reaction  Thermochemistry 

In  Tables  lA-IC  are  displayed  our  QCISD(T)  calculated  (and,  where  known, 
the  experimental)  values  for  the  relative  energies  of  the  reactant  M'*'(ns^;  ^S)  + 
H2,  excited-  state  reactant  M+fnsnp;  ^’*P)  +  H2,  and  product  MH'*’  +  H  and  HMH'*' 
species.  In  all  cases,  the  energies  are  derived  from  electronic  energies;  no 
zero-  point  corrections  are  included. 

The  lowest  excited  and  ^P  sutes  of  M***  are  listed  because  they  give  rise 
to  excited  ^*^82.  and  ^.iBi  stales  of  C2v  MH2''’  which,  in  turn,  affect  the 

ground-  state  reaction  dynamics  via  second-  order  Jahn-  Teller  coupling^  ^  to 
or  intersections  with  the  ground  ^Ai  state  as  described  later  in  this  paper.  It  is 
essential  that  our  calculations  place  these  excited  states  reasonably  accurately 
if  our  inferences  about  the  ground-  state  dynamics  are  to  be  reliable. 

An  important  point  to  note  about  these  data  is  that  the  experimental 
thre.<holds  for  producing  MH(D)'*’  *  H(D)  do  not  correlate  with  the 
thermodynamic  energy  differences 


AExhermo.  =  E(MH(D)+)  +  E(H(D))  -  E(M+)  -  E(H2  (D2  or  HD)) 


which  appear  in  the  fifth  rows  of  Tables  I.  Nor  do  these  thresholds  apree  with 
‘'e  locations  of  the  "barriers"  on  the  potential  energy  surfaces  shown  in 
Figs.  1  and  listed  in  the  seventh  rows  of  Tables  I.  These  facts  make  it  clear  that  a 
"dynamical"  rather  than  energetic  constraint  must  be  operative  in 
determining  the  experimental  thresholds  which  exceed  by  from  ca.  0.4  eV  (for 
BD'*’  formation)  to  as  much  as  5  eV  (for  GaD'*’  formation)  the  thermodynamic 
energy  requirements.  It  is  for  this  reason  that  we  must  now  turn  our  attention 
to  the  dynamics  of  the  M*  +  H-H  collision. 


IV.  Dynamics 


A.  Experimental  Conditons  and  Their  Implications 


1.  Initial  Conditions 

The  guided  ion  beam  experiments  of  Aimentrout  and  co-workers^ 
involve  collisions  in  which  the  H2  (D2  or  HD)  and  M***  reagents'  internal 
(vibrational,  rotational,  and  electronic)  degrees  of  freedom  usually  exist  in  or 
close  to  thermal  equilibrium  near  room  temperature.  Therefore,  nuclear 
motions  along  these  degrees  of  freedom  are  restricted,  in  the  early  stag- 
ion-  molecule  collision,  to  narrow  ranges  approximately  characterized  by  the 
corresponding  classical  turning  points.  As  a  result,  the  most  important  areas  of 
the  potential  energy  surface  in  the  entrance  channel  region  are  those  for 
which  such  internal  modes  do  not  deviate  greatly  from  their  most  probable 
values. 

In  contrast,  the  relative  kinetic  energy  betweeen  the  M'*'  ion  and  its  H2 
(or  D2  or  HD)  collision  partner  is  very  large  in  comparison  with  thermal 
energies.  This  collision  energy,  and  its  associated  momentum,  has  components 
along  three  directions:  (i)  the  M-to-H  (or  D)  axis  (rMH)*  (>>)  ibe  other  M-to-H  (or 
D)  axis  (rMH')>  (iiO  t^e  out-of-molecular  plane  angular  coordinate 
Explicitly,  the  classical  collisional  kinetic  energy  in  an  M***  +  A-B  encounter  is 


^collision 


1  mM(mA-*-mB)  dY  2  (ny^+n^ 

2  mA-t-me-t-mM  ^dt^  ~  ^“mA+mB+mM’ 


where  the  kinetic  energy  of  the  M***  ion  as  measured  (and  experimentally 
controlled)  in  a  lab-  fixed  coordinate  system  is 


Ttab  - 


MM  dY  2 
2  Mt^ 


Here.  mM.  inA>  and  mg  represent  the  masses  of  the  three  particles,  and  Y  the 
separation  of  M  to  the  center  of  mass  of  the  A-B  pair.  The  collisional  kinetic 
energy  can  be  decomposed  into  components  describing  motion  of  M'*’  along  the 
fMA  and  r^B  axes  as  follows: 


Tma  =  Ttab 


BA  - 
mA+mM 


=  T  collision 


OtA 

mA+mg 


.BA  Am  A.  2 

=  2  ^  dt  ^ 


TmB  = 


TLab 


mg 


mg+mM 


=  ^collision 


mg 

mA-t-mg 


drM£  2 

=  2  ^  dt  ^ 


where  the  second  equalities  are  only  approximate  because  they  assume  mM  » 
mA  mg. 

2.  Role  of  Kinetic  Energy  Along  Collision  Degrees  of  Freedom 

The  above  decomposition  of  the  collision  energy  has  been  used^  to 
rationalize  the  occurrence  of  different  energy  thresholds  for  production  of 
MH'*'  and  MD'*'  in  M'*’  HD  reactions,  the  idea  being  that  there  is  more  energy 
along  the  rMD  axis  (2/3  Tcollision)  tiian  along  the  rMH  axis  (1/3  Tcollision).  so 
MD'*’  can  be  formed  at  lower  total  collision  energies.  However,  this  model 
predicts  that  the  difference  in  thresholds  for  MH'*’  and  MD'*'  should  be  large 
(with  the  MD'*'  threshold  occuring  at  one  half  the  collision  energy  of  MH'*).  and 
that  for  M"*  +  H2  or  M"*  +  D2,  where  there  are  1/2  Tcollision  and  2/4 
Tcollision  along  the  rj^H  and  rMD  axes,  the  thresholds  should  be  even  different 
than  in  the  HD  case.  These  quantitative  deuils  are  not  observed^  in  the 
experimental  data,  although  there  are  significant  differences  (much  more 
than  zero-  point  energies  can  account  for)  in  the  HD,  H2,  and  D2  thresholds. 

The  primary  difficulty  with  using  this  fractional  collision  energy  Tm  B 
concept  is  that  it  ignores  how  the  potential  energy  V  depends  on  the  two 
independent  rMH  and  rMD  axes.  If  V  were  a  strong  function  of  one  of  these 


COO',  ifuates  (c.g.,  of  the  distance  from  M'*'  to  the  nearest  H  or  D  center),  and 
depended  weakly,  if  at  all,  on  the  other  distance,  then  this  decomposition  of 
Tcollision  would  make  more  sense.  When  the  kinetic  energy  along  the 
"important"  M-to-H  (or  D)  coordinate  was  adequate  to  overcome  any  barrier  in 
V  along  this  same  coordinate,  reaction  could  take  place.  However,  for  the 
reactions  at  hand,  V  depends  on  the  two  rviH  and  rMH'  (or  rMD)  distances  in  a 
symmetrical  fashion;  that  is,  the  electronic  energy  remains  the  same  if  rM  H 
and  tMH'  are  interchanged.  Moreover,  V  is  a  strong  function  of  both  distances, 
at  least  within  the  entrance  channel  where  the  collision-  to-  internal  energy 
transfer  is  initiated.  As  a  result,  kinetic  energy  along  both  rMH  and  n^H' 
required  to  access  regions  of  the  potential  where  reaction  can  occur;  neither 
^MH  Qor  Tmh'  alone  is  adequate.  It  is  for  these  reasons  that  consideration  of  the 
kinetic  energy  alone  does  not  adequately  explain  the  isotope  effects  on 
thresholds. 

Nevertheless,  the  different  masses  of  the  H  and  D  isotopes  do,  in  fact, 
have  important  effects  on  the  thresholds  for  MD'*'  and  MH'*'  formation,  but  not 
because  of  the  reasons  outlined  above.  The  Hydrogenic  masses  mA  and  ms,  as 
well  as  the  H-D,  H2,  or  D2  reduced  mass  p,  appear  in  the  kinetic  energy, 
approximately  (see  Sec.  IV.  G  for  more  detail)  as: 

T  -  HA .  BA a.  dr  2 

‘  =  2  ^  dt  ^  2  ^  dt  ^  ^  2  W  • 


As  discussed  later  in  this  ptqter,  so-  called  mass-  weighted  coordinates  r'  =  'Vp  r, 
r'MD  =  Vmo  rMD  and  ^mh  =  V m h  rMH  can  be  introduced  after  which  the  total 
energy  H  =  T  V  is  expressed  as: 


1  ,<ir'M A.2  ,  1  .dr'MB.2  , 
2^  dt  ^  "^2^  dt  ^ 


2  '^^'■’ma.  r’MB.  r). 


In  this  form,  isotopic  differences  disappear  from  T  and  appear  only  in  the 
different  dependence  of  V  on  r'MA<  ^MB>  and  r'.  Although  V  depends  on  rMA  and 
rMB  is  a  symmetrical  manner,  its  dependence  on  r'MA  and  r'MB  loay  be 
asymmetric  and  reflects  the  A  and  B  masses.  It  is  these  mass  dependences  that 
produce  isotope  effects  in  the  local  normal-  mode  frequencies  derived  from 
such  a  Hamiltonian,  and  it  is  these  mass  effects  that  we  think  more  correctly 
explains  the  isotope  effects  on  reaction  thresholds. 


B.  Entrance-  Channel  Reaction  Dynamics 


The  potential  euergy  function  along  the  relative-  motion  degrees  of 
freedom  is  slowly  varying  as  the  collision  begins  (i.e.,  at  large  R  and  small  r). 
As  the  collision  progresses,  these  three  degrees  of  freedom  evolve  in  a  manner 
that  produces  significant  forces  (i.e..  changes  in  potential)  along  rMH  i^MH 
Keeping  in  mind  that  essentially  all  of  the  initial  momentum  is  directed  along 
these  "soft  modes",  and  recalling  that  restoring  forces  strive  to  preserve  C2v 
symmetry,  we  direct  attention  to  flux  moving  with  high  initial  energy  and 
velocity  along  both  rMH  ^nd  rMH‘  liule  energy  along  the  r-  axis  (because  of 
the  low  vibrational  energy  of  the  H2  reagents). 

C.  Entering  the  Region  of  Strong  Interaction 

As  flux  progresses  up  the  entrance  channel  to  higher  potential  energy, 
the  kinetic  energy  and  momenta  along  imh  and  rMH*  decrease,  but  lack  of 
coupling  between  the  R-  and  r-  directions  (i.e.,  (d^E/dRdr)  «  0  as  illustrated 
clearly  in  Figs.  1)  makes  energy  (and  momentum)  transfer  from  the  relative- 
motion  modes  to  the  transverse  r-  dominated  mode  very  ineffective. 

However,  as  flux  moves  to  even  smaller  R-  values,  a  region  of  space  is 
reached  where  energy  transfer  can  occur.  This  region  is  characterized  not 
only  by  curvature  (i.e.,  existence  of  off-  diagonal  d^E/dRdr  terms)  on  the 
potential  energy  surfaces  as  shown  in  Figs.  1  but  also  by  near  degeneracies  in 
the  eigenvalues  of  the  mass-  weighted  Hessian  (MWH)  matrix  (see  below) 
evaluated  at  such  geometries. 

These  statements  now  need  to  be  justified  by  introducing  and  using  the 
MWH  matrix  as  a  device  for  analyzing  the  dynamical  resonances  that  permit 
energy  transfer  and  subsequent  chemical  reaction  to  occur. 

D.  The  Hessian  as  a  Local  Approximation  to  the  Potential  Energy 

The  Hessian  matrix,  evaluated  at  a  geometry  in  the  region  of  strong 
interaction  (denoted  {x|(<’})  and  expressed  in  terms  of  the  3N  Cartesian 
coordinates  (xk)  of  the  N  atoms  is 


H'k.m  =  (3^E/dxk3xni)xjjO  . 


The  gradient  vector 

F'k  =  OE/dxk)x^o 


evaluated  at  this  same  point  gives  the  slope  of  the  energy  along  the  Cartesian 
directions  xk-  Of  course,  the  values  of  this  matrix  and  vector  depend  strongly  on 
where  these  derivatives  are  evaluated;  at  a  point  (xk°)  in  the  strong  interaction 
region,  {F’k)  has  large  components  along  the  inter-  fragment  coordinates. 

These  constructs  allow  the  potential  energy  surface  V(xk)  to  be 
approximated  to  the  point  {xk°}  as  a  Taylor  series: 

V(xk)  =  V({xkO))  +  Ik  F’k  6xk  +  1/2  Ikjn  H’k.m  5xk  8xni. 

where  5xk  means  the  deviation  of  xk  from  the  value  xk°. 

E.  The  Kinetic  Energy  in  Mass  Weighted  Coordinates 

Of  course,  the  kinetic  energy  T  can  also  be  written  in  terms  of  the  3N 
Cartesian  displacement  coordinates  {5xk}-  However,  if  so-  called  mass  weighted 
coordinates 

yk  =  Xk 

are  introduced,  where  mk  is  the  mass  of  the  nucleus  to  which  the  coordinate  xk 
pertains,  the  kinetic  energy  can  be  written  as  a  sum 

of  3N  terms  each  of  which  has  the  same  (unit)  mass  factor.  By  so  treating  the 
kinetic  energy  in  a  manner  that  assigns  equal  mass  to  all  3N  degrees  of 


freedom,  the  potential  energy  function  alone  governs  the  natural  frequencies 
of  motion  of  the  system. 

F.  The  Mass  Weighted  Hessian 

In  terms  of  these  {yic}  coordinates,  the  local  quadratic  approximation  to 
the  potential  energy  is  given  by: 

V(yk)  =  2k  Fk  8yk  +  1/2  2k.m  Hk.m  8yk8ym  . 

where  Fk  is  the  gradient 

Fk  =  OE/ayk)x^o 

=  (mk)*l/2  (dE/axk)x,^o 

of  the  electronic  energy  along  the  yk  coordintate,  Hk,m  is  the  matrix  of  second 
derivatives  with  respect  to  the  yk  variables: 

Hkjn  «  O^E/aykdym)x,jO 
=  H'k,m  (“k  mm)*  1/2 

and  5yk  is  the  displacement  along  the  yk  coordinate  from  the  point  at  which  the 
derivative  is  evaluated.  The  matrix  {Hk.m)  is  called  the  mass-  weighted  Hessian 
matrix  (MWH),  and  (Fk)  is  the  gradient  vector  in  mass-  weighted  coordinates. 
Notice  that  Hv  m  has  units  of  sec'2  because  5yk  has  units  of  gml/2  cm; 
therefore,  the  eigenvalues  of  (Hk.m)  introduced  in  the  next  section  have  units 
of  sec'2,  or  frequency  squared. 

G.  The  Classical  Equations  of  Motion 

A  classical  Hamiltonian 

H.T  +  V»Zkj  +  Ik  Fk  5yk  +  1/2  Ikm  Hk.m  SykSym 


treatment  can  be  used  to  describe  the  (local)  motion  of  the  3N  degrees  of 
freedom.  The  Newton  equations  of  motion  then  read 


d^Syic 

dt2 


-Fk-I^m  Hk  8ym  ■ 


The  linear-  plus-  quadratic  form  of  the  potential  is  a  reasonable 
representation  of  the  potential  along  internal  degrees  of  freedom  of  the 
fragments  (e.g..  the  r-  coordinate)  because  these  modes  undergo  small- 
amplitude  motions  about  their  equilibrium  positions.  However,  this  is  an 
unreasonable  global  representation  for  the  potential  along  inter-  fragment 
degrees  of  freedom.  The  latter  coordinates  are  not  bounded  by  the  potential  at 
large  R,  whereas  the  quadratic  terms  above,  if  (Hk.m)  is  positive  definite,  would 
constrain  the  inter-  fragment  coordinates  at  large  separations.  The  model 
potential  thus  produces  harmonic  oscillatory  motion  even  along  the  inter¬ 
fragment  degrees  of  freedom,  which,  of  course,  really  undergo  collisions  with  a 
single  close  encounter  between  the  fragments.  Nevertheless,  as  shown  below, 
the  description  of  inter-  fragment  degrees  of  freedom  (i.e.,  rMH  and  rMH') 
provided  by  these  equations  of  motion  is  useful  in  analyzing  the  dynamics  local 
to  the  points  {xk°}  of  strong  interaction  and  for  the  short  duration  of  the 
collision. 


1.  The  MWH  Eigenmode  Basis 

To  make  further  progress,  we  now  introduce,  for  reasons  that  will  soon 
become  clear,  the  (unitary)  matrix  uk,j  that  diagonalizes  the  MWH  matrix  {Hk.ml 

Hjyn  Umj  *  Uk,j» 


and  we  denote  the  non-  zero  eigenvalues  by  (j  =  1,  2 .  3N-5  or  3N-6).  The 

MWH  matrix  will  also  have  S  or  6  eigenvalues  and  corresponding  eigenvectors 
belonging  to  the  translation  and  rotation  of  the  entire  Mnn***  specie.  We  denote 
these  vectors,  which  can  be  formed  independent  of  diagonalizing  the  MWH 
matrix,  by  {tk.j;  j  =  1,  2 . 5  or  6), 


If  the  Newton  equations  are  now  multiplied  by  ukj  and  summed  over  k, 
we  obtain  equations  of  motion 


d^SAj 
dt2  ■ 


-fj  -  (0j2  8Aj 


for  the  components  5Aj  of  5y)c  along  the  normalized  eigenmodes  of  the  MWH; 
5Aj  =  2k  Uk.j  5yk. 

Here,  fj  is  the  projection  of  the  {Fk )  force  vector  along  the  eigenmode  of  the 
MWH: 


fj  =  2k  Uk.j  Fk  . 


2.  The  MWH  Model  Dynamical  System  for  Bimolecular  Dynamics 


The  equations  derived  above 


-fj  -  (0j2  8Aj 


specify  the  time  evolution  of  a  model  dynamical  system  containing  3N-S  or 
3N  -6  modes  that  undergo  sinusoidal  motions  (if  the  Q)j2  are  positive) 


8Aj(t)=  5Aj(eq)  +  Aj(t^)  cos(a)jt) 


at  frequencies  (dj  (sec*^)  about  equilibrium  positions 


5Aj(eq)  = 


Here,  SAj(t=0)  is  the  amplitude  of  motion  along  the  normal  mode,  which  is.  in 
turn  related  to  the  total  energy  Ej  contained  in  that  mode: 

Ej  =  ^  l^^^2  4.  i  (o2j  l5Aj  -  8Aj(eq)|2 


=  2  IAj(t=0)|2  [  sin2((Dji)  +  cos2(o)jt)  ]. 


so 


Let  us  now  examine  how  this  model  dynamics  relates  to  the  M***  -<•  A-B  collision 
dynamics  under  study. 

a.  Internal  Modes 

In  the  example  at  hand,  forces  along  internal  modes  of  the  ion  or  its 
collision  partner  (i.e.,  the  r-  dominant  H-H  vibration)  are  small  in  the  early 
stages  of  the  collision  because  of  the  small  excursions  experienced  by  these 
degrees  of  freedom,  so  oscillatory  motion  does  indeed  take  place  about  the 
equilibrium  position.  Also,  the  energy  content  of  these  modes  is  small,  so  the 
corresponding  amplitudes  Aj(tsO)  will  be  small  and  can  be  estimated  as 

where  kT  is  the  thermal  energy.  For  these  modes,  the  picture  provided  by  the 
MWH  model  is  appropriate  and  easily  understood. 

b.  Relative-  Motion  Modes 

In  contrast,  the  MWH  picture  of  the  inter-fragment  motions  (i.e.,  the 
modes  arising  from  rMH  >iid  rMH')  requires  further  examination.  At  points  {xk°} 
on  the  potential  energy  surface  where  strong  coupling  between  the  rMH  or 
rMH'  ^d  internal  (r)  modes  are  likely,  the  forces  fj  along  inter-  fragment 
coordinates  will  be  large  and  repulsive  (see  Figs.  1).  The  curvature  of  the 
potential  surface  along  these  directions,  as  reflected  in  the  corresponding 
eigenvalues  of  the  MWH.  will  be  positive  (see  later)  and  substantial. 


The  MWH  dynamical  model  treats  these  inter-  fragment  degrees  of 
freedom  as  also  undergoing  harmonic  motion,  but  about  a  minimum  that  is  far 

■  fj  fj2 

removed  (by  an  amount  — 7  )  from  the  point  {xk®}  and  which  lies  ~^r  in 

energy  below  its  value  at  {xk°}.  Clearly,  this  description  of  the  inter-  fragment 
motion  is  not  globally  correct  because  the  true  collisional  dynamics  invloves  a 
single  encounter  between  the  fragments,  not  a  sinusoidal  series  of  such 
enounters.  Nevertheless,  if  used  only  for  the  brief  time  interval  during  which 
strong  mode  coupling  is  realized,  this  does  gives  a  useful  local  model  of  the  true 
dynamics  because: 


(i)  it  adequately  describes  the  potential  surface  (i.e.,  the  forces  and  local 
natural  frequencies  of  motion)  near  points  {xj<)}  where  mode  coupling  is 
strongest,  and 


(ii)  it  includes  the  correct  relative  kinetic  energies  along  all  modes. 


For  these  reasons,  the  approximate  MWH  Newton  equations  can  be  used  to 
obtain  the  time  evolution  of  the  system  for  the  (brieO  duration  of  the  collision 
during  wnich  the  ion  resides  in  this  repulsive  region  of  the  potential 
surface  and  during  which  energy  transfer  is  possible. 

H.  Avoided  Crossings  of  MWH  Eigenvalues 

At  geometries  where  a  (local)  relative*  motion  MWH  eigenvalue  o),  and 
an  internal*  mode  eigenvalue  osint  undergo  an  avoided  crossing,  there  is 
enhanced  probability  of  energy  transfer  from  the  collision  coordinate  to  the 
mode  associated  with  mint.  In  such  cases,  one  can  think  of  the  dynamics  of  two 
coupled  oscillators  having  frequencies  0)3  and  (Oint>  whose  coordinates  obey 

^^2  *  ~  ®  ®^int 


where  £1  (with  units  sec"^)  denotes  the  coupling  between  the  two  coordinates. 
In  the  absence  of  coupling,  these  two  coordinates  would  undergo  simple 
sinusoidal  motions  about  their  own  equilibrium  positions  and  at  their  own 
frequencies. 

However,  as  shown  in  many  elementary  classical  mechanics  texts,  ^  ^ 
when  coupling  is  present,  the  time  evolution  involves  two  new  characteristic 
frequencies  (i)±  .  In  the  limit  where  (0$  =  cDin^  m  coq  (i.e.,  when  the  two  natural 

frequencies  would  cross  if  coupling  were  absent),  the  two  new  frequencies  are 
given  by: 

(i>±  =  >/o)o^  ±  O, 


which  reduces  to 

“±  =  “0*2^ 

if  im  «  (00^.  as  it  is  for  the  cases  considered  below. 

If  energy  E  is  initially  deposited  into  the  collision  mode  s.  the  the  two 
modes  will  evolve  in  time  as: 


6As(t)  -  5As(eq)  =  -yj  ^  )  cos(o)ot) 

8Aint(t)  -  5Aint(eq)  = '\j  ^  ^  sin(<Dot) 


For  short  times  t,  the  SAs(t)  mode  contains  much  energy  and  its  coordinate 
oscillates  at  frequency  (dq-  As  time  evolves,  the  5Aint(t)  mode  gains  amplitude, 
and  once 

Qt 

2<no  ”2 

this  mode  has  acquired  all  of  the  amplitude  (and  hence  energy)  that  the  5As(t) 

nmp  ^ 

mode  originally  had.  One  thus  says  that  in  a  time  interval  t  =  "q  the  energy 
transfer  takes  place;  alternatively,  the  rate  of  energy  transfer  is: 


18 


rate 


a 

Jtwo  ■ 


This  result  would  be  most  relevant  if  the  coupling  £2  were  operative  as 
detailed  above  throughout  the  entire  sinusoidal  motions  of  the  two  oscillators. 
However,  to  model  the  situation  at  hand,  it  is  more  proper  to  allow  £2  to  act  only 
for  the  narrow  range  of  inter-  fragment  distances  AR  where  the  two  modes 
undergo  their  avoided  crossing.  A  modification  of  the  above  rate  expression 
that  allows  £2  to  act  only  for  that  fraction  f  of  an  oscillation  (of  sin((oot))  that 
the  collision  resides  within  AR  is  given  as  follows: 

£^  .  £2  AR  £2AR 

rate  = -  f  = - (on  = - 

K(oo  K(oo  V  nv 


Here  ~  is  the  residence  time  of  the  trajectory  with  speed  v  in  the  range  AR, 

and  (00  is  the  inverse  of  the  time  it  takes  to  make  one  oscillation.  Of  course,  the 
speed  v  can  be  expressed  in  terms  of  the  energy  E  in  the  s  mode,  and  the 
potential  V  at  the  geometry  where  the  avoided  crossing  occurs. 

A  substantial  body  of  experience  in  the  classical  dynamics  of  multi-  mode 
systems  has  shown  that  when  the  (local)  natural  frequencies  of  two  degrees  of 
freedom  become  nearly  equal  (  <0$=  <‘Hnt  ■  (ao  )•  energy  transfer  between  these 
modes  is  most  likely.  Within  a  quantum  dynamics  treatment,  energy  transfer  is 
facile  when  two  modes  have  equal  or  nearly  equal  energy  spacines.  The 
classical  and  quantum  points  of  view  are  easily  seen  to  be  consistent  when,  as 
here,  a  local  quadratic  treatment  (which  incorporates  the  true  local  forces  and 
curvatures)  is  used  for  the  potential.  In  such  a  case,  the  resultant  harmonic 
frequencies  a>t  <»int  8>ve  both  the  natural  frequencies  of  the 
corresponding  periodic  motions  and  the  frequency  spacines  between 
neighboring  quantum  states  that  differ  by  unit  quantum  number.  Thus  the 
resonance  condition  discussed  above  can  be  viewed  either  as  near  equality 
between  two  natural  periodic  oscillation  times  or  as  near  equality  between  two 
quantum-  state  energy  spacings. 

It  is  also  known  that  movement  through  regions  of  such  near 

ar 

degeneracy  must  have  a  "conuct”  or  "residence"  time  ~  long  enough  to 
permit  the  coupling  between  the  two  modes  that  unergo  the  avoided  crossing  to 


effect  a  transition.  If  movement  through  this  region  is  extremely  fast,  energy 
transfer  is  unlikely.  In  the  following  Section,  such  avoided  crossings  are  used 
to  explore  under  what  conditions  such  energy  transfer  can  readily  occur. 

V.  Findings  and  Relation  to  Experiments 

A.  Avoided  Crossings 

We  show  the  eigenvalues  of  the  locally  calculated  MWH  for  M"*"  =  B'*’.  Al‘*‘. 
and  Ga'*'  in  Figs.  (4A-4C),  and.  in  each  case,  results  for  all  three  isotopes  (H,  HD, 
and  DD)  are  shown.  In  Fig.  4D  the  eigenvalues  of  the  MWH  are  shown  for  the 
collinear  approach  path  for  comparison.  In  all  cases,  the  distance  r  between 
the  two  hydrogenic  centers  was  held  fixed  at  the  equilbrium  value  in  H2.  0.7SS 
A.  This  was  done  because  the  geometries  that  play  critical  roles  in  determining 
where  energy  transfer  occurs  are  not  those  in  which  all  nuclear  coordinates 
are  "relaxed”,  but  those  that  would  be  realized  during  high-  energy  ion- 
molecule  collisions  such  as  those  taking  place  in  the  guided  ion  beam 
experiments.  At  least  in  the  early  stage  of  such  collisions,  before  energy 
transfer  has  taken  place,  the  H-H  (or  H-D  or  D-D)  distance  deviates  only  slightly 
from  0.755  A. 

In  all  of  Figs.  4,  the  relative-  motion  eigenvalues  are  very  small  at  large- 
R  where  the  forces  between  M***  and  H2  (or  D2  or  HD)  are  quite  weak,  and  the 
internal-  mode  eigenvalue  is  large.  As  R  decreases,  the  former  eigenvalues 
increase  because  the  inter-  fragment  forces  increase,  and  eventually  one  or 
more  avoided  crossings  (or  actual  crossing  for  the  H2  and  D2  cases  in  which  the 
asymmetric  stretching  mode  is  uncoupled  by  symmetry  from  the  two  ai  modes) 
take  place. 

The  energy  transfer  ideas  reviewed  above  imply  that  facile  energy  (and 
momentum)  transfer  from  the  (soft)  rMH  imh'  collision  eigenmodes  into  the 
r-  dominated  internal  mode  can  occur  near  an  avoided  crossing.  i£  a  collision 
has  enough  kinetic  energy  to  access  these  avoided  crossing  regions.  From  Fig. 
4D.  which  pertains  to  the  collinear  geometry  case,  we  note  that  avoided 
crossings  do  not  occur  at  all,  at  least  within  the  energy  range  studied.  This 
combines  with  the  bending  mode's  geometric  instability  of  the  linear 
structures  to  further  emphasize  the  importance  of  near-  C2v  geometries 
relative  to  near-  collinear  geometries. 


B.  Relation  to  Reaction  Thresholds 


For  all  of  the  species  considered  here,  as  shown  in  Figs.  4A-  4C,  the 
avoided  crossings  occur  at  geometries  where  the  potential  energy  is 
considerably  in  excess  of  either  the  thermodynamic  threshold  or  the  "barrier" 
on  the  C2v  potential  surface.  In  Table  III,  the  inter-  fragment  distances  (R)  at 
which  the  avoided  crossings  occur  (i.e.,  where  the  splitting  between 
interacting  MWH  eigenvalues  are  smallest)  are  listed  as  are  the  potential 
energies  at  these  geometries.  The  experimental  thresholds  for  formation  of 
MH'*’  and  MD'*',  where  known,  are  also  listed. 

It  should  be  noted  that  the  interactions  among  modes  that  gives  rise  to 
the  avoided  crossings  do  not  exist  only  at  the  R-  values  listed  in  Table  III.  Such 
interactions  are  present  over  a  significant  range  of  inter-  fragment  distances, 
and  certainly  develop  significant  strength  somewhat  before  reaching  the  R- 
values  listed.  For  this  reason,  we  specify  lower  bounds  to  the  critical 
interaction  distances  when  we  quote  geometries  where  the  MWH  eigenvalues 
come  closest.  Moreover,  because  the  potential  energy  surfaces  are  quite  "steep" 
and  repulsive  in  these  regions,  the  energies  derived  at  our  quoted  R-  values 
represent  upper  bounds  to  the  minimum  energies  needed  to  effect  reaction. 

Having  made  these  qualifying  remarks,  the  model  dynamics  provided  by 
the  MWH  eigenmode  analysis  explains  the  reaction  thresholds  in  terms  of  the 
avoided  crossings.  When  the  kinetic  energy  of  collision  Tcollision  large 
enough  to  access  geometries  where  the  MWH  eigenvalues’  undergo  avoided 
crossings,  energy  transfer  to  the  internal  mode  (r)  induces  reaction.  As  inter¬ 
fragment  collisional  kinetic  energy  is  lost,  energy  is  deposited  into  the 
internal  mode,  thereby  causing  the  H-H  (D-D  or  H-D)  bond  to  lengthen  and  to 
eventually  rupture.  The  geometries  at  which  these  avoided  crossings  occur  are 
typified  by  strong  repulsive  forces  along  both  rMH  and  rMH*  (or  rMD)  axes. 
Therefore  collisions  that  access  these  regions  must  have  high  kinetic  energies 
along  both  of  these  axes.  For  this  reason,  it  is  the  total  kinetic  energy,  not  Tm  a 
or  Tmb>  tbat  is  the  key  collision  energy  parameter. 

The  data  summarized  in  Table  III  clearly  show,  for  example,  that 
thresholds  for  B***,  Al*^,  and  Ga'*’  reacting  with  D2  should  occur  ca.  1.3  eV,  2.5  eV. 
and  3.3  eV  above  their  respective  endothermicities.  Although  our  predicted 
thresholds  display  trends  much  like  the  experimental  findings  (see  Table  III). 


the  energies  where  the  avoided  crossings  are  strongest  tend  to  systematically 
exceed  the  experimental  thresholds  by  ca.  1  eV  (see  comments  above  about 
upper  bounds).  This  is  likely  a  result  of  the  steeply  repulsive  nature  of  the 
potentials  (e.g..  the  energies  drop  by  more  than  1  eV  over  a  0.05  A  range  of  R  m 
these  regions  for  all  three  species)  at  such  geometries  and  the  fact  that 
significant  mode  coupling  develops  at  longer  R-  values  than  where  the  avoided 
crossing  is  strongest.  In  addition,  the  thermal  motions  of  the  H2,  D2.  or  HD  are 
non-  zero  and  tend  to  make  the  apparent  experimental  thresholds  lower  than 
the  true  thresholds. 

C.  Isotope  Effects  for  HD 

In  the  HD  cases,  the  two  relative-  motion  modes  have  different  natural 
frequencies;  the  mode  dominated  by  rvfH  motion  has  higher  frequency  than 
that  dominated  by  rMD-  a  result,  the  former  mode  undergoes  an  avoided 
crossing  with  the  internal  (r-  dominated)  mode  at  larger  R,  and  hence  at  lower 
energy.  Energy  that  is  thus  transferred  from  the  rMH  motion  to  the  internal 
mode  decreases  the  relative  velocity  along  rMH.  but  not  (as  much)  along  rMD- 
The  differential  velocity  that  thus  develops  between  rMD  aod  rMH  causes  the  M- 
to-D  distance  to  shorten  more  rapidly  than  the  M-to-H  distance,  while  the  H-to- 
D  distance  is  growing  (since  energy  is  being  put  into  this  mode  to  break  the  H- 
D  bond).  As  these  movements  propagate  in  time,  MD'*'  is  formed  and  H  is 
eliminated. 

The  important  point  is  that  the  lower-  energy  avoided  crossing  involves 
coupling  energy  out  of  the  rMH  mode  and  production  of  MD"**  +  H.  Likewise,  the 
higher-  energy  avoided  crossing,  which  involves  the  rMD*  dominated  mode 
coupling  to  the  H-D  motion,  produces  MH'*'  +  D.  The  difference  in  thresholds  for 
MD'*’  and  MH'*'  formation  is  explained  by  differences  in  the  energies  at  which 
fbe  rMH  tmd  rMD  avoided  crossings  occur. 

D.  Coupling  Strengths 

The  avoided  crossing  graphs  also  provide  information  about  the  strength 
of  coupling  between  the  relative-  motion  and  internal  modes.  When  the 
eigenvalues  oi^±  that  "avoid"  one  another  are  viewed  as  solutions  of  a  2  x  2 

matrix  eigenvalue  problem,  the  difference  ((i>^.t.  -  (o^.)  between  them  can  be 


related  to  the  off-  diagonal  element  of  the  matrix  (which  we  denote  n  and 
which  has  units  of  sec'^): 

o  -  (0^.) 

“  "  2 

In  Table  III,  we  also  report  these  coupling  strengths  Vn  (in  cm'^  units)  for  all 
of  the  cases  considered  here.  Thus  far.  we  have  not  made  a  direct  connection 
between  these  coupling  strengths  and  experimental  findings.  Clearly,  they 
relate  to  the  magnitude  of  energy  flow  between  the  relative-  motion  and 
internal  modes  near  the  avoided  crossing,  but  their  magnitudes  do  not  seem  to 
correlate  with  ion  yield  or  branching  ratio  (e.g..  ii  is  larger  for  coupling  to  the 
fMH  mode  than  to  the  rMD  mode,  although  the  yield  (of  MD***)  from  the  former  is 
smaller  than  for  the  latter).  Of  course,  the  observed  ion  yields  are  not  direct 
measures  of  the  initial  rate  of  formation  of  MD'*'  or  MH***  because  nascent  ions 
may  undergo  decomposition  before  being  detected,  and  the  fraction  that 
decompose  depends  on  the  collision  energy  E. 

E.  Long  Interaction  Times  are  Required 

Not  only  must  a  "trajectory"  access  the  avoided  crossing  geometry,  it 
must  spend  enough  time  there  to  permit  the  couplings  to  effect  energy 
transfer.  The  time  (t)  spent  in  this  region  can  be  estimated  in  terms  of  the 
initial  collision  energy  Ecoii.  the  range  of  R-  values  over  which  the  coupling 
takes  place  AR,  the  electronic  potential  energy  near  the  avoided  crossing  V^ross 
as  well  as  the  reduced  mass  of  the  M'*’  -i-  A-B  pair 

gAa..”M 

^  fflAB  +  oiM 
AR 

V2(Ecoll  -  VcrossVlt 


The  time  needed  to  effect  energy  transfer  is  related  to  the  strength  of 
coupling  between  the  two  modes  undergoing  the  avoided  crossing.  As  shown 
earlier,  this  coupling  (Q  in  sec'^)  can  be  extracted  from  the  avoided  crossing 


graphs  as  one  half  the  "splitting"  between  the  two  eigenvalues  at  their  closest 

.  „  (0)2+  -  0)2.) 
approach;  Q  = - - -  . 


So,  if 


X  Vn  =  —  -  (O^low) 

^  V2(Ecoll  *  Vcj.oss)/M 

energy  transfer  can  be  facile.  This  implies  that  collsions  with  incident  kinetic 
energies  slightly  in  excess  of  V^j^oss  t>e  most  effective  in  transferring 
energy  into  the  t-  dominated  degrees  of  freedom,  and  that  collisions  with  much 
higher  kinetic  energy  should  be  less  effective.  For  the  cases  considered  here, 
as  shown  in  Table  III,  I/V2  V((o^high  *  w^low)  ranges  from  ca.  952  cm'^  to 
1953  cm*^  (i.e.,  corresponding  to  frequencies  of  3-6  xl0^3  sec*^)  and  AR  ranges 
from  0.05  A  to  ca.  0.2  A.  Therefore,  one  expects  collisions  passing  through  the 
avoided  crossing  region  at  ca.  10^  to  10^  cm  sec'^  or  slower  to  be  quite  effective. 
This  means  that  collisions  with  kinetic  energy  along  the  collision  mode  to 
which  r  is  coupled  much  in  excess  of  the  potential  at  the  avoided  crossing  will 
be  ineffective. 

F.  Fate  of  Collisions  that  Result  in  Energy  Transfer 

Those  collisions  that  access  geometries  where  energy  transfer  from  a 
relative-  motion  coordinate  to  an  internal  mode  can  occur  have  a  chance  to 
evolve  to  produce  MH'*’(or  MD'*')  product  ions.  In  doing  so,  the  H-H  (D-D  or  H-D) 
bond  breaks,  a  new  M-H  (or  M-D)  bond  is  formed,  and  an  H  (or  D)  is  eliminated. 

In  the  picture  provided  by  the  MWH  eigenmode  model,  once  enough 
energy  and  momentum  are  transferred  to  the  r-  coordinate,  flux  can  evolve 
toward  larger  r-  values.  Such  flux  will  move  toward  the  "barrier"  regions  of  the 
potential  energy  surfaces  shown  in  Figs.  1,  although  the  total  energy  exceeds 
the  barrier  energy  (of  3.2,  4.5,  and  4.6  eV  for  B*,  AI'*’.  and  Ga"*",  respectively)  by 
more  than  1  eV  in  all  cases. 

However,  as  flux  so  evolves,  our  analysis  of  the  three  potential  surfaces 
in  Figs.  1  shows  that  a  region  on  the  ^  A 1  potential  surface  is  reached  within 
which  either  (i)  the  asymmetric  stretch  motion  of  b2  symmetry  becomes 


unstable  (i.e..  develops  a  negative  MWH  eigenvalue)  due  to  second-  order  Jahn- 
Teller  coupling  with  the  nearby  *82  excited  state  or  (ii)  the  ^82  excited  state 
intersects  and  passes  below  the  ^Ai  surace.  In  either  case,  flux  can  move,  with 
no  restoring  forces,  away  from  C2v  symmetry.  It  is  this  step  that  permits  the 
asymmetric  rupture  of  the  M species  to  produce  the  observed  MH'*’  (or  MD'*’) 

+  H  (orD).  In  Tables  II  are  shown  the  geometries  at  which  the  ^82  state  has  its 
own  minimum  because  it  is  near  such  geometries  that  the  second-  order  Jahn- 
Teller  couplings  or  surface  intersections  are  most  likely.  Also  shown  in  Tables 
II  are  the  eigenvalues  of  the  MWH  at  the  "barrier"  geometry;  in  all  cases,  one 
notes  an  imaginary  frequency  for  the  b2  mode,  which  reflects  the  geometrical 

instability  of  these  regions  to  asymmetric  distortion. 

8ecause  the  regions  of  avoided  crossings  of  MWH  eigenvalues  occur  high 
above  the  "barrier"  regions,  the  MH'*'  (or  MD'*')  product  ions  are  likely  to  be 
formed  with  a  large  amount  of  internal  (vibration/rotation)  and  translational 
energy.  8ecause  the  M-H'*'  bond  strengths  are  rather  weak  (48,  18,  and  IS 
kcal/mol  for  HH*^,  AIH*^,  and  GaH'*'.  respectively),  such  internal  energy  can 
cause  the  nascent  MH*^  species  to  fragment  before  reaching  the  experiment's 
detector.  Hence,  fragmentation  of  the  product  ions  can  contribute  to  the 
unusually  small  cross-  sections^  found  for  these  reactions,  although  another 
cause  is  likely  to  severe  "steric"  requirments  imposed  by  reaching  the  region 
of  strong  mode  coupling  and  the  inefficient  relative-  motion  to  interna! - 
motion  energy  flow. 

VI.  Summary 

Energies  at  which  the  local  natural  frequencies  corresponding  to  inter¬ 
fragment  and  to  internal  motions  (obtained  as  eigenvalues  of  the  MWH  matrix) 
undergo  avoided  crossings  are  related  to  kinetic  energy  thresholds  in  the  ion- 
molecule  reactions  M'*’  +  H2  -♦  MH"*"  +  H,  for  M"*"  =  B"*",  Al'*’,  and  Ga'*’  and  deuterium 
substituted  analogs.  It  is  the  total  kinetic  energy  of  collision  TcoUision  • 
components  Tma  Tmb  along  the  two  M-to-H  (or  D)  axes,  that  determines 
whether  a  collision  can  access  the  regions  of  the  potential  surface  where 
avoided  crossings  occur.  Although  systemmatic  differences  exist  between  the 
(lower  bound)  apparent  experimental  thresholds  and  our  (upper  bound) 
predictions,  the  trends  seem  to  be  in  agreement.  Moreover,  the  fact  that 


thresholds  exceed  thermodynamic  requirements  is  reproduced  by  our  model,  as 
is  the  propensity  to  produce  MD'*'  at  lower  collision  energy  than  MH"*". 

The  primary  assumption  in  making  correlations  between  reaction 
thresholds  and  avoided  crossings  of  the  MWH  eigenvalues  is  that  energy 
transfer  in  such  mode-  coupling  collisions  is  the  rate  determining  step  in 
forming  MH'*'  products.  Such  a  model  was  introduced  because  the 
experimentally  observed  reaction  thresholds  exceed  by  0.4  to  ca.  5  eV  the 
thermodynamic  energy  requirements  or  computed  barrier  heights  of  these 
reactions  and  because  the  measured  cross-  sections  are  very  small.  This  is.  of 
course,  not  true  for  all  ion-  molecule  reactions,  but  is  for  the  "impulsive" 
reactions  considered  here. 

The  fact  that  the  avoided  crossings  occur  high  above  the  thermodynamic 
thresholds  leads  to  large  internal  energies  in  .li,.  (or  MD'*')  product  ions 

and  to  likely  subsequent  dissociation.  This  is  one  of  the  likely  causes  for  the 
measured  cross-  sections  for  MH'*'  (or  MD***)  formation  being  small  (smallest  for 
Ga"*"  and  largest  for  B"*’). 

Preference  to  form  MD"*"  at  lower  collision  energies  than  MH'*’  when  HD 
reacts  with  M***  is  consistent  with  the  avoided-  crossing  frequency-  resonance 

picture  introduced  here.  The  higher  frequency  M-H  mode  (which  leads  to  MD"*" 
products)  couples  to  the  high  frequency  internal  motion  (H-D)  mode  at  larger 

R-  values  (and  hence  lower  energy)  than  the  lower  frequency  M-D  mode. 

Although  the  MWH  matrix  is  used  as  a  tool  in  this  analysis,  the  model  put 
forth  here  is  not  equivalent  to  a  reaction  path  Hamiltonian^  dynamics  model, 
which  also  employs  the  MWH.  The  latter  introduces  a  reaction  path  connecting 
a  transition  state  (i.e.,  a  first-  order  saddle  point  on  the  energy  surface)  to  the 
reagent  geometry.  The  critical  geometries  of  our  approach  (those  where 
avoided  crossings  of  MWH  eigenvalues  occur)  have  energies  much  in  excess  of 
the  nearest  first-  order  saddle  points,  and  thus  do  not  lie  on  a  reaction  path.  Our 
geometries  relate  more  to  paths  that  would  be  realized  in  high-  energy  ion- 
molecule  collisions  in  which  the  reagents  have  little  internal  energy. 
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Figure  Captions 

Fig  lA.  C2v  Symmetry  Contour  Plot  of  the  (^Ai)  Ground  State  Energy  of  B"*"  +  H2. 
The  R  (distance  of  B"*"  to  the  center  of  H-H)  and  r  (H-H  distance)  axes  are  in  A, 
and  the  contours  are  spaced  by  10.0  kcal/mol.  The  symbol  X  is  used  to  denote 
the  location  of  the  "barrier",  and  Y  is  used  to  denote  the  region  of  strong  mode 
mixing  (see  text). 

Fig  IB.  C2v  Symmetry  Contour  Plot  of  the  (*Ai)  Ground  State  Energy  of  Al'*'  +  H2. 

The  R  (distance  of  Al"*"  to  the  center  of  H-H)  and  r  (H-H  distance)  axes  are  in  A, 
and  the  contours  are  spaced  by  10.6  kcal/mol.  The  symbol  X  is  used  to  denote 
the  location  of  the  "barrier",  and  Y  is  used  to  denote  the  region  of  strong  mode 
mixing  (see  text). 


Fig  1C.  C2v  Symmetry  Contour  Plot  of  the  ('Ai)  Ground  State  Energy  of  Ga"*"  +  H2. 
The  R  (distance  of  Ga*  to  the  center  of  H-H)  and  r  (H-H  distance)  axes  are  in  A, 
and  the  contours  are  spaced  by  10.4  kcal/mol.  The  symbol  X  is  used  to  denote 
the  location  of  the  "barrier”,  and  Y  is  used  to  denote  the  region  of  strong  mode 
mixing  (see  text). 


Fig.  2  Coordinate  System  Used  to  Label  C2v  Geometries. 

Fig.  3.  Contour  Plot  of  the  (^Ai)  Ground  State  Energy  of  B'*’  +  H2  for  collinear 
geometries.  The  r(MH)  (distance  of  B+  to  the  nearest  H  atom)  and  r  (H-H 
distance)  axes  are  in  A,  and  the  contours  are  spaced  by  6.3  kcal/mol. 

Fig.  4A  Avoided  Crossings  of  Eigenvalues  of  the  Mass-  Weighted  Hessian  Matrix 
for  B'*’  +  H2.  D2>  and  HD.  The  horizontal  axis  is  R  (A)  and  the  vertical  axis  is  o) 
(cm* I).  For  large  R,  the  hipest  frequency  mode  is  the  HH,  DD,  or  HD 
stretching  vibration,  and  the  lower  two  are  the  relative-  motion  modes. 

Fig.  4B  Avoided  Crossings  of  Eigenvalues  of  the  Mass-  Weighted  Hessian  Matrix 

for  Al'*'  -f  H2.  D2.  and  HD.  The  horizontal  axis  is  R  (A)  and  the  vertical  axis  is  o) 

(cm*^).  For  large  R,  the  highest  frequency  mode  is  the  HH,  DD,  or  HD 
stretching  vibration,  and  the  lower  two  are  the  relative-  motion  modes. 

Fig.  4C  Avoided  Crossings  of  Eigenvalues  of  the  Mass-  Weighted  Hessian  Matrix 

for  Ga'*'  H2.  D2,  and  HD.  The  horizontal  axis  is  R  (A)  and  the  vertical  axis  is  u 

(cm*^).  For  large  R,  the  highest  frequency  mode  is  the  HH,  DD,  or  HD 
stretching  vibration,  and  the  lower  two  are  the  relative-  motion  modes. 

Fig.  4D.  Plot  of  Eigenvalues  of  the  Mass-  Weighted  Hessian  Matrix  for  B'*'  H2. 

D2.  and  HD  in  collinear  geometries  with  the  H-H  stretch,  interfragment,  and 
bending  vibrations  labeled.  The  horizonul  axis  is  R  (A)  and  the  vertical  axis  is 
(D  (cm*^). 


Table  lA.  Electronic  States  Energies  (kcal/mol)  Measured  with  Respect  to 

B+(^S)  +  H2('2:g‘*'). 


Species* 

This  work*’ 

Experiment* 

B+(^S)  + 

0.0 

0.0 

B+(^S)  +  2H(^S) 

109 

110 

B+(^P)  +  H2(^2:/) 

107 

107 

B+(^P)  + 

215 

210 

BH+  (2l)  +  H 

AEThermo.  =  *5 1 

61 

HBH+(^I/) 

-  60 

B+  -H2  barrier 

73 

a.  C  E.  Moore.  Tables  of  Arnmtc  Energy  Levels,  and  K.P.  Huber,  G.  Her2berg, 
Molecular  Spectra  and  Molecular  Structure.  IV.  Constants  of  Diatomic 
Molecules  (Van  Nostrand  Reinhold,  New  York.  1979). 

b.  Based  on  QCISD(T)  dau  except  for  Ip  sute  where  projected  fourth-  order 
M0Uer-  Plesset  (PMP4)  perturbation  theory  was  used  due  to  difficulties  in  the 
QCISD(T)  convergence. 


Table  IB.  Electronic  States  Energies  (kcal/mol)  Measured  with  Respect  to 

Ai+(^S)  +  H2(‘2g'^). 


This  work^ 

Experiment^ 

A1+(^S)  +  H2(*2/) 

0.0 

0.0 

A1+(4)  +  2H(^S) 

109 

110 

A1+(^P)  +  H2(^2/) 

105 

107 

A1+(^P)  +  H2(^Z/) 

180 

171 

A1H+  (21)  +  H 

AE-rhermo.  =  91 

HA1H+(^Z/) 

12 

A1'*'—H2  barrier 

104 

a.  C.  E.  Moore.  Tables  of  Atomic  Energy  Levels,  and  K.P.  Huber,  G.  Herzberg, 
Molecular  Spectra  and  Molecular  Structure.  IV.  Constants  of  Diatomic 
Molecules  (Van  Nostrand  Reinhoid.  New  Yoric.  1979). 

b.  Based  on  QCISD(T)  data  except  for  state  where  projected  fourth'  order 
M0ller-  Plesset  (PNff4)  perturbation  theory  was  used  due  to  difficulties  in  the 
QCISD(T)  convergence. 


Table  1C.  Electronic  States  Energies  (kcal/mol)  Measured  with  Respect  to 

Ga+{^S)  +  H2(‘5:g'*'). 


Species* 

This  work*’ 

Experiment* 

Ga+(‘S)  +  H2(‘I/) 

0.0 

0.0 

Ga+l^S)  +  2H(^S) 

109 

110 

Ga+(^P)  +  H2(^I/) 

123 

137 

Ga+(^P)  +  H2(^2:/) 

196 

202 

GaH+  (21)  +  H 

^EThermo.  » 

HGaH+(^le‘^) 

20 

Ga+"  H2  barrier 

105 

a.  C.  E.  Moore,  Tables  of  Atomic  Energy  Levels,  and  K.P.  Huber,  G.  Herzberg, 
Molecular  Spectra  and  Molecular  Structure.  IV.  Constants  of  Diatomic 
Molecules  (Van  Nostrand  Reinhold,  New  York,  1979). 

b.  Based  on  QCISD(T)  data  except  for  state  where  projected  fourth-  order 
M0ller-  Plesset  (PMP4)  perturbation  theory  was  used  due  to  difEculties  in  the 
QCISD(T)  convergence. 


Table  IIA.  Total  and  Relative  Energies.  Geometries,  and  Vibrational  Frequencies 
for  Species  Relating  to  the  B'*'  +  H2  BH'*'  +  H.  HBH'*’  Reactions 


Species 

Electronic 

Energies 

(Hartrees) 

Optimized 
Internuclear 
Distances  (A) 

Vibrationai 
Frequencies^/ 
Zero  Point 
Energies 
(cm-1 ) 

Relative 

Energies 

(kcal/mol)3 

B+(1S)  +  H2 

-25.446250 

-25.468830 

— 

4224/21  12 

0.0 

0.0 

BH+(2i)  +  H 

-25.351373 

-25.372139 

rBH  =  1-199 

2582/ 

1291 

60 

61 

HBH+(1z+g) 

-25.520364 

-25.564074 

r.  2.374 

2594  (ai), 
2880  (b2). 

932  (bend)/ 
3669 

-47 

-60 

B+---H2  barrier 

-25.322627 

-25.352085 

wm 

78 

73 

BH2+  (IB2) 
minimum 

-25.328460 

-25.399644 

r=1.674 

R-0.996 

1026  (ai), 
2173(b2), 
2083(ai)/ 
2641 

74 

43 

a.  In  all  cases,  the  energies  are  given  relative  to  the  B***  +  H2  reactants.  These 
are  electronic  energies,  and  thus  do  not  include  zero-point  corrections.  In 
each  case,  and  for  the  column  giving  total  energies  in  Hartrees,  the  first 
number  is  based  on  our  CAS-MCSCF  calculations,  and  the  second  is  based  on  our 
QCISD(T)  darn. 

b.  These  local  harmonic  frequencies  were  obtained  from  the  analytical  second 
derivatives  of  the  MCSCF  energy  at  the  MCSCF  geometries. 


Table  IIB.  Total  and  Relative  Energies.  Geometries,  and  Vibrational  Frequencies 
for  Species  Relating  to  the  Al'*’  +  H2  ==>  AlH"*"  +  H.  HA1H+  Reactions 


Species 

Electronic 

Energies 

(Hartrees) 

Optimized 
Internuclear 
Distances  (A) 

Vibrational 
Frequencies*^/ 
Zero  Point 
Energies 
(cm'1 ) 

Relative 

Energies 

(kcal/mol)3 

AI+(1S)  +  H2 

-242.856705 

-242.879646 

r-  0.755 

4224/2112 

0.0 

0.0 

AIH+(22)  +  H 

-242.717076 

-242.735419 

rAIH  -  1-658 

1424/ 

712 

88 

91 

HAIH+(lZ+g) 

-242.804625 

-242.860474 

r-  3.103 

1940(31), 

2055  (b2). 
513(bend)/ 

251 1 

33 

12 

AI+---H2  barrier 

-242.692731 

-242.713641 

r  -  1.852 

R  «  1.587 

2362  i  (ai) 

1942  i  (b2) 

996  (ai) 

103 

104 

AIH2+  (IB2) 
minimum 

-242.685868 

-242.740915 

r-1.729 

R-1.429 

855  (ai), 
1307(b2), 
1637(ai)/ 
1900 

107 

87 

a.  In  all  cases,  the  energies  are  given  relative  to  the  Al*  -f  H2  reactants.  These 
are  electronic  energies,  and  thus  do  not  include  zero-point  corrections.  In 
each  case,  and  for  the  column  giving  total  energies  in  Hartrees,  the  Hrst 
number  is  based  on  our  CAS-MCSCF  calculations,  and  the  second  is  based  on  our 
QCISD(T)  dau. 

b.  These  local  harmonic  frequencies  were  obtained  from  the  analytical  second 
derivatives  of  the  MCSCF  energy  at  the  MCSCF  geometries. 


33 


Table  IIC.  Total  and  Relative  Energies.  Geometries,  and  Vibrational  Frequencies 
for  Species  Relating  to  the  Ga'*’  +  H2  =>  GaH"*"  +  H,  HGaH'*'  Reactions 


Species 

Electronic 

Energies*^ 

(Hartrees) 

Optimized 
Internuclear 
Distances  (A) 

Vibrational 
Frequencies*^/ 
Zero  Point 
Energies 
(cm-'* ) 

Relative 

Energies 

(kcal/mol)3 

Ga+(iS)  +  H2 

-258.1 19452 
-1924.206695 

r-  0.757 

4224/21  12 

0.0 

0.0 

GaH+(2Z)  +  H 

-257.967714 

-1924.056193 

^GaH  -  t.747 
roaH  -  1.65 

902/451 

83 

94 

HGaH+(iI+g) 

-258.069058 

■1924.174511 

r  -  3.096 

2003(ai), 

2139(b2), 

628(bend)/ 

2699 

32 

20 

Ga'*'—H2  barrier 

r  -  2.0 

R  -  1.75 

not  available*^ 

105 

GaH2‘*'  (^82) 
minimum 

•257.909068 

-1923.991427 

r.1.886 

R.1.390 

132 

135 

a.  In  all  cases,  the  energies  are  given  relative  to  the  Ga***  +  H2  reactants.  These 
are  electronic  energies,  and  thus  do  not  include  zero-point  corrections.  In 
each  case,  and  for  the  column  giving  total  energies  in  Hartrees,  the  first 
number  is  based  on  our  CAS-MCSCF  calculations,  and  the  second  is  based  on  our 
QCISD(T)  dau. 

b.  These  local  harmonic  frequencies  were  obtained  from  the  analytical  second 
derivatives  of  the  MCSCF  energy  at  the  MCSCF  geometries. 

c.  The  MCSCF  calculations  used  a  pseudo-  potential,  but  the  QCISD(T)  data 
involve  all  electrons. 

d.  The  finite  difference  routines  used  in  GAMESS  were  not  able  to  produce  a 
reliable  Hessian  matrix  in  this  case. 


Table  III.  Geometry,  Energy,  and  Coupling  Strength  in  the  Region  of  Avoided 
Crossing  of  Mass  Weighted  Hessian  Eigenvalues, 
and  Experimental  Reaction  Thresholds 


Species 

R  at 

Crossing  (A) 

AEThermo 

/eVb 

E  at  Crossing 
(kcal/mol;  eV) 

Coupling 

Strength 

Vn  (cm*l) 

Experimental® 
Thresholds  (eV) 
to  form  (MA'*’) 

B+  +  HH 

>1.05 

2.6 

<  89:  3.9 

1953 

3.3  ±0.1 

B+  +  DD 

>1.05 

<  89;  3.9 

1302 

3.3  ±0.1 

B+  +  HD 

>  1.00  (rMD) 

>  1.05  (rviH) 

<  107: 4.6 

<  89:  3.9 

1121 

1502 

4.0  ±  0.2  (BH+) 

3.0  ±  0.2  (BD+) 

A1++  HH 

>  1.22 

3.9 

<  148:  6.4 

1352 

6.6  ±0.2 

A1+  +  DD 

>  1.22 

<  148;  6.4 

952 

6.6  ±0.1 

A1+  +  HD 

>1.16  (rMD) 

>  1.22  (rMH) 

<  177:  7.7 

<  148:  6.4 

1070 

1121 

6.7  ±0.1  (AlH-*-) 

4.7  ±0.1  (A1D+) 

Ga+  +  HH 

>  1.21 

4.1 

<  170:  7.4 

1881 

not  available 

Ga+  +  DD 

>1.21 

<  170: 7.4 

1323 

8.5  ±  0.5  (GaD+) 

Ga+  +  HD 

>  1.15  (rMD) 

>  1.25  (rMH) 

<  206;  9.0 

<  145;  6.3 

1120 

1483 

not  available 

a.  Ref.  1  as  well  as  P.  B.  Armentrout  (private  communication). 

b.  Not  zero-  point  corrected,  so  independent  of  isotopic  masses. 
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